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ABSTRACT Heteroresistance refers to phenotypic heterogeneity of microbial clonal populations under antibiotic stress, and it has 
been thought to be an allocation of a subset of "resistant" cells for surviving in higher concentrations of antibiotic. The assump- 
tion fits the so-called bet-hedging strategy, where a bacterial population "hedges" its "bet" on different phenotypes to be selected 
by unpredicted environment stresses. To test this hypothesis, we constructed a heteroresistance model by introducing a 
b/a CTX _ M _ 14 gene (coding for a cephalosporin hydrolase) into a sensitive Escherichia coli strain. We confirmed heteroresistance in 
this clone and that a subset of the cells expressed more hydrolase and formed more colonies in the presence of ceftriaxone (ex- 
hibited stronger "resistance"). However, subsequent single-cell-level investigation by using a microfluidic device showed that a 
subset of cells with a distinguishable phenotype of slowed growth and intensified hydrolase expression emerged, and they were 
not positively selected but increased their proportion in the population with ascending antibiotic concentrations. Therefore, 
heteroresistance — the gradually decreased colony-forming capability in the presence of antibiotic — was a result of a decreased 
growth rate rather than of selection for resistant cells. Using a mock strain without the resistance gene, we further demonstrated 
the existence of two nested growth-centric feedback loops that control the expression of the hydrolase and maximize population 
growth in various antibiotic concentrations. In conclusion, phenotypic heterogeneity is a population-based strategy beneficial 
for bacterial survival and propagation through task allocation and interphenotypic collaboration, and the growth rate provides a 
critical control for the expression of stress-related genes and an essential mechanism in responding to environmental stresses. 

IMPORTANCE Heteroresistance is essentially phenotypic heterogeneity, where a population-based strategy is thought to be at 
work, being assumed to be variable cell-to-cell resistance to be selected under antibiotic stress. Exact mechanisms of heteroresis- 
tance and its roles in adaptation to antibiotic stress have yet to be fully understood at the molecular and single-cell levels. In our 
study, we have not been able to detect any apparent subset of "resistant" cells selected by antibiotics; on the contrary, cell popu- 
lations differentiate into phenotypic subsets with variable growth statuses and hydrolase expression. The growth rate appears to 
be sensitive to stress intensity and plays a key role in controlling hydrolase expression at both the bulk population and single-cell 
levels. We have shown here, for the first time, that phenotypic heterogeneity can be beneficial to a growing bacterial population 
through task allocation and interphenotypic collaboration other than partitioning cells into different categories of selective ad- 
vantage. 
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The increasing antibiotic resistance of bacterial pathogens is of eroresistance: phenotypic heterogeneous resistance to an antibi- 
great concern worldwide. The exploration of possible mecha- otic treatment within a clonal population (1-3). In laboratory 
nisms of such resistance has led to intensive investigations of het- practice, such resistant subpopulations emerge as resistant colo- 
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FIG 1 A plasmid-bearing ESBL gene leads to heteroresistance. (a) PAPs of four clinical isolates. Four isolates were analyzed, including two E. coli (Eco) strains 
(carrying the CTX-M-14 and TEM-10 ESBL genes) and two K. pneumoniae (Kpn) strains (producing CTX-M-65 and SHV-2). Ceftriaxone was the antibiotic 
used for all strains, (b) Construction of the plasmid expressing CTX-M-14. (c) Fluorescent image of cells from overnight culture of the resistant strain R4. 
Growing cells show significant differences in CTX-M-14-GFP expression. The inset shows a histogram of GFP intensity for the same cell population. 



nies grown inside the inhibitory halos in Etest and lead to a grad- 
ual decrease in the colony count on agar plates as the antibiotic 
concentration increases in a population analysis profile (PAP) as- 
say (4). 

Heteroresistance is observed in a wide range of microbes (5-7) 
and often in cases of treatment failures (8, 9); it has been thought 
to be a prerequisite of resistant subpopulations that lead to full 
resistance due to their selective advantages upon antibiotic stress 
(10-12). However, this speculation remains controversial because 
the "resistant" colonies do not always give rise to resistant proge- 
nies or converge to a homotypic fully resistant phenotype after 
antibiotic selection (7). In addition, the genetic basis of heterore- 
sistance remains obscure, although attempts to discover it have 
been made (13, 14) and a single resistant gene, mecA, has been 
reported to be responsible for heteroresistance in Staphylococcus 
aureus (15). 

Although heteroresistance generally refers to a subset of phe- 
notypically resistant cells from antibiotic-sensitive clinical iso- 
lates, phenotypic heterogeneity in terms of antibiotic resistance 
has often been observed in resistant strains. In fact, a sensitive or 
resistant strain is recognized according to artificial criteria, and 
the heteroresistance of these strains may share common mecha- 
nisms that vary only with the resistance mode. Here, our study was 
focused on a particular case where a hydrolase was used to con- 
struct a simple cellular model: a plasmid containing a gene encod- 
ing an extended-spectrum beta-lactamase (ESBL) (a beta-lactam 
hydrolase), CTX-M-14, was introduced into a commonly used 
laboratory antibiotic-sensitive Escherichia coli strain, DH5a. We 
subsequently investigated population-based and single-cell-based 
performance of the resistant strain under broadly varying antibi- 
otic stresses (concentrations). The hydrolase gene appears to be 



tightly coupled with growth rate control; antibiotic levels, resis- 
tance gene expression, and cell growth form feedback loops exhib- 
iting growth plasticity and adaptation to the changing antibiotic 
stress. 

RESULTS AND DISCUSSION 

Host heteroresistance reflects heterogeneous expression of hy- 
drolase genes. To see if ESBL genes are related to heteroresistance 
and using PAP, we analyzed four ESBL-bearing clinical isolates 
from Klebsiella pneumoniae and Escherichia coli, all of which ex- 
hibit heteroresistance in the presence of ceftriaxone (Fig. la). The 
behavior of these strains is in accordance with a previous report 
that a bla OXA-58-like gene-encoded hydrolase was demonstrated 
to always lead to phenotypic heterogeneous resistance (PHR) 
(heteroresistance) to carbapenems in Acinetobacter baumannii 
(1). 

To carry the investigation further, we first constructed a mo- 
lecular model for heteroresistance by introducing a plasmid 
(pUA66) carrying an ESBL gene, encoding CTX-M-14, into E. coli 
DH5a. We also inserted the CTX-M-14 gene promoter upstream 
of the green fluorescence protein (GFP) open reading frame 
(ORF) (16) (Fig. lb) so that the fluorescence intensity of GFP 
directly reflected the expression of the CTX-M-14 gene. The rea- 
son we chose the CTX-M-14 gene for the study is based on its 
widely known prevalence (17, 18) and the fact that there is no trans 
regulator or cis elements controlling its expression when antibiot- 
ics are in action (19). Second, we evaluated heteroresistance and 
CTX-M-14 gene expression in the resistant host strain, namely, 
E. coli R4. The PAP curve of R4 for ceftriaxone appears similar to 
that of the original clinical E. coli isolate from which the CTX- 
M-14 gene was cloned, with a slight shift toward the left (Fig. la). 
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FIG 2 High-intensity group (HIG) and low-intensity group (LIG) cells show 
differences in resistance to antibiotics, (a) Sorting gates of HIG and LIG cells, 
(b) Etest of LIG (left) and HIG (right) cells, (c) Copy numbers of the plasmid 
and CTX-M-14 transcripts in HIG and LIG cells, (d) Dynamic GFP intensity 
histograms of HIG and LIG cells after sorting. Cells of both subgroups were 
cultured in antibiotic-free medium and sampled for FCM assay every 30 min 
for 4 h immediately after sorting. The red dotted line indicates the peak posi- 
tion of the whole population under antibiotic- free conditions, and the blue 
dotted lines indicate peak positions of the sorted LIG and HIG subpopulations. 



The GFP intensity of the R4 overnight culture had a roughly log- 
arithmic Gaussian distribution in flow cytometry (FCM), ranging 
from 10 2 to 10 5 arbitrary units (AU), with a peak at 1,238 AU 
(Fig. lc). Consequently, the CTX-M-14 gene itself was heteroge- 
neously expressed and sufficiently led to detectable heteroresis- 
tance at the same time. 

Hydrolase gene expression in subpopulations correlates pos- 
itively with resistance and adaption to variable antibiotic stress. 
To further confirm if hydrolase gene expression correlates with 
resistance, we flow sorted (fluorescence-activated cell sorting 
[FACS] ) the top 7% (named the high-intensity group [HIG] ) and 
the lowest 10% (named the low- intensity group [LIG]) bacterial 
populations based on their GFP intensity from an overnight cul- 
ture grown in drug-free medium (Fig. 2a). We used Etest to mea- 
sure the ceftriaxone MICs of the LIG and HIG, which are 4 ju,g/ml 
and >256 /xg/ml, respectively (Fig. 2b); this result indicates that 
higher GFP intensity or more CTX-M-14 hydrolase production 
corresponds to stronger resistance. To see if plasmid copy number 
variation contributes to the expression of the hydrolase, we quan- 
tified copy numbers of both the plasmid and the CTX-M-14 tran- 
script in the two GFP intensity groups, using real-time PCR and a 
housekeeping gene, gapA, as a reference; the corresponding fold 
changes in the copy numbers of the plasmid and the hydrolase 
transcript for HIG over LIG are 4.89X and 7.67X, respectively 
(Fig. 2c). The different fold changes reflect variations in plasmid 
copy number and transcription activity between the two sub- 
groups. 

Further observation uncovered that the hydrolase gene exhib- 
ited quantized expression when we cultured the HIG and LIG cells 
in antibiotic-free medium immediately after sorting and investi- 
gated the dynamic distribution of the expression. The LIG- 



associated peak at 732 AU eventually switched to the position of 
the antibiotic- free peak at 1,238 AU in -30 to 60 min. The initial 
HIG-associated peak at 23,158 AU faded away in ~4 h and was 
replaced by the antibiotic-free peak at 1 ,238 AU, which emerged at 
-60 min (Fig. 2d). In both cases, the peaks did not shift continu- 
ously but were replaced by peaks characteristic of the antibiotic- 
free condition. 

We further examined the antibiotic effect with hydrolase ex- 
pression using exponential-phase cells incubated in a series of 
double dilutions of ceftriaxone for 3 h. Along with the increased 
ceftriaxone dosage, cell proliferation was suppressed (Fig. 3a). 
Meanwhile, the average florescence intensities and MIC in Etest 
for the entire population shared a common trend, ascending grad- 
ually but with a constant shoulder at the dosage interval of 8 to 
128 ju.g/ml (Fig. 3b). The result further confirms the correlation 
between CTX-M-14 expression and resistance. 

The distribution of GFP intensities appeared to be fixed at 
several values, which further confirms the discontinuity of CTX- 
M-14 hydrolase gene expression. When the antibiotic concentra- 
tion varied from 2 to 128 ju-g/ml, fluorescence peaked at 4,580 AU. 
When the antibiotic concentration reached a higher level, 
256 /^g/ml and above, several peaks appeared, with a major one at 
13,473 AU (Fig. 3c). This observation suggests the existence of 
complex subpopulations, where their CTX-M-14 hydrolase gene 
expression is heterogeneous and their proportions adapt to anti- 
biotic dosages. 

In addition, our microscopic observation confirms that cellu- 
lar fluorescence intensity increases with ceftriaxone concentra- 
tions and cells grow in a filament-like way in response to higher 
concentrations (Fig. 3d) . However, flow cytometry can detect only 
the total fluorescent intensity of the bacterial particles, which are 
often elongated to be severalfold longer than the untreated cells. 
Although forward scatter (FSC) reflects particle size in flow cy- 
tometry, it is not exactly proportional to particle size, and we 
cannot normalize fluorescence to a single cell unit by simply di- 
viding by the FSC value. 

Increased CTX-M-14 hydrolase expression does not provide 
cells with a selective advantage but is associated with a reduced 
growth rate. According to the exclusion rule (12), heteroresis- 
tance can be a result of heterogeneous expression of antibiotic 
hydrolase genes and provides a heterogeneous selective advantage; 
when treated with antibiotics, the cells with higher resistance are 
selected and form resistant colonies in inhibitory halos. To test 
this rule at the single-cell level, we fabricated a PDMS (polydim- 
ethylsiloxane) microfluidic device that allows tracing a single bac- 
terial cell and recording cellular activity under different antibiotic 
stresses. Fluorescence (GFP) intensity can be recorded at each site 
of the tiled cells instead of the total value of the entire particle in 
flow cytometry. We started the experiment with R4 culture in 
exponential phase; 50 randomly selected cells at each ceftriaxone 
concentration were monitored for 2 h to observe their colony 
formation, GFP intensity, and growth rate. 

We found that almost all cells, even dim cells — expressing 
CTX-M-14 at a lower level — had excellent viability and displayed 
heterogeneity in growth rate and GFP intensity at various antibi- 
otic concentrations below 2,048 ju,g/ml. We did, however, observe 
cell death at an extreme ceftriaxone concentration, 8,196 )LLg/ml, 
when some cells died as determined by fluorescence fade-out or 
sudden burst after cell wall bulging. Surprisingly, at any antibiotic 
concentration, brighter cells (higher GFP/CTX-M-14 expression) 
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FIG 3 Effects of antibiotic on the GFP intensity of the R4 strain, (a) Growth curves of the R4 strain cultured in medium containing serially diluted ceftriaxone, 
(b) Average fluorescence intensity (left y axis; a.u., arbitrary units) and MIC measured with Etest (right y axis) at the time point of 3 h of ceftriaxone treatment 
as shown in panel a. (c) Histograms of GFP intensity after 3 h of ceftriaxone treatment as shown in panel a. They axis in these histograms indicates the percentage 
of the maximum cell count (% Max), and the color lines indicate the ceftriaxone concentrations (u.g/ml) used. Dotted lines labeled with "L" (low), "M" 
(moderate), or "H" (high) indicate the positions of the major peaks, which shift rightward as the ceftriaxone concentration increases, (d) Fluorescence images of 
cells after 2 h of ceftriaxone treatment as shown in panel a. 



did not show a significantly higher survival rate than dim ones 
(Fig. 4a). Although brighter cells often grow slower, and slow- 
growing cells are statistically more resistant to antibiotics 
(Fig. 4b), as previously reported (20), we did not observe a direct 
relationship between cellular GFP intensity and the survival rate: 
even the brightest cells often faded out, a sign of death. We further 
constructed a mock plasmid without the CTX-M- 14 ORF, but the 
GFP ORF remained downstream of the CTX-M- 14 gene pro- 
moter. Due to the loss of the resistance gene, the sensitive strain 
(named S4) could sustain only ceftriaxone concentrations lower 
than 128 ng/ml, when cells death became apparent. As expected, 
the S4 cell survival rate at a ceftriaxone concentration of 128 ng/ml 
showed no relationship with GFP intensity but negatively corre- 
lated with the growth rate (Fig. 4c and d). 

When cells were treated with a series of ascending ceftriaxone 
dosages which the R4 strain can sustain, the mean and variance of 
the GFP intensity of the cell population both increased (Fig. 5a), 
which is in line with our observations in the flow cytometry ex- 
periment. Another major effect of antibiotic treatment was a grad- 
ual reduction of the growth rate (Fig. 5a). This observation is in 



accordance with the growth curves of the bulk cultures subjected 
to the same series of antibiotic treatments. We have now identified 
a unique process in R4 when it responds to antibiotic stress. First, 
the majority of the cells are fast growing, expressing a low level of 
GFP, i.e., a low-expression phenotype. Second, when subjected to 
different dosages of ceftriaxone, a subset of cells switch their orig- 
inal phenotype to high expression and slow growing, while others 
maintain their original low-expression and fast-growth pheno- 
type. Third, the percentage of phenotype-switching cells increases 
with an elevated antibiotic dosage (Fig. 5b). 

We also noticed a phenotype switch happening during the split 
of a mother cell into two daughter cells (CI and C2; Fig. 5e); one 
(CI) of the daughter cells retained its original growth rate, with a 
relatively slow accumulation of cellular GFP, while C2 acquired 
another trend: high expression and slow growth. There seems to 
be a control on total activity of cellular growth and hydrolase 
expression. 

Let us reexamine the concept of heteroresistance, which origi- 
nated with colony-based assays (Etest and PAP). In fact, colony 
formation is not an indicator of individual living cells inoculated 
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FIG 4 Survival rates of resistant (R4) and control sensitive (S4) strains at bactericidal ceftriaxone concentrations, (a) Correlation of death rates of R4 cells with 
their initial GFP intensities after 2 h of ceftriaxone treatment at the minimal bactericidal concentration (8,192 fig/ml). Cells were binned according to their initial 
GFP intensities, (b) Correlation of death rates of R4 cells with their initial growth rates after the same treatment as for panel a; cells were binned according to their 
initial growth rates, (c and d) Graphs corresponding to panels a and b for S4 treated at its minimal bactericidal concentration of ceftriaxone (128 ng/ml). Death 
count (black squares) and survival rate (blue squares) were plotted together. 



on the agar plate, but colony-forming potential depends on dif- 
ferences between the growth rate and the death rate in a probabi- 
listic manner. According to our results at the single cell level, an- 
tibiotic concentrations used in such assays do not select for cells 
with "resistant" phenotypes but lead to a general effect of in- 
creased CTX-M- 14 expression and decreased growth rate through 
phenotype switching toward heterogeneity. The gradual reduc- 
tion of the growth rate decreases the possibility of colony forma- 
tion and leads to a decreased colony count as the antibiotic dosage 
increases. Therefore, a clear picture emerges as we take the growth 
rate into account. On agar plates with antibiotic, some cells switch 
to the high-expression and slow-growing phenotype and hydro- 
lyze antibiotic in its surrounding space, creating favorable condi- 
tions for their fast-growing progenies (although the progenies 
may have the low-expression phonotype) to form colonies. Since 
it is not a survival game but a growth game where a clonal popu- 
lation is living through antibiotic stress, we hypothesize that a 
population-based growth strategy is intrinsic to bacterial species 
through phenotype diversification and cooperation among prog- 
enies, and thus the subpopulations sacrifice some of their energy 
and other resources, such as replication, to CTX-M- 14 gene ex- 
pression, achieving maximal population survival and propaga- 
tion. 

Phenotypic heterogeneity of clonal populations has been 
found to be a common phenomenon for a long time (21-23). 
Heterogeneity is deemed to offer progenies with different selective 
advantages in case of unpredictable environmental stresses, such 
as the persister phenotype under extreme antibiotic stress (24, 25). 
The population has to hedge its bet on different phenotypes and 



hope that at least one of the phenotypes may survive when the 
others perish. However, in a harsh environment, it may not be 
easy for a strain to develop a single phenotype that bears all the 
characteristics needed to survive that stress; an effective alterna- 
tive strategy is of essence for the population to allocate different 
tasks to heterogeneous subpopulations, each with its own capaci- 
ty/phenotype based on collaborative strategies and opportunity 
gains among the subpopulations. In our study, we identified two 
opposite phenotypes, high expression and slow growth being one 
and low expression and fast growth the other; neither secured life 
independently, but corporately the whole population managed to 
obtain its maximal biomass under stresses that are often variable 
and certainly unpredictable. 

Adaptation to variable antibiotic dosages is achieved by 
stringent control of growth-centric feedback loops. Since CTX- 
M- 14 expressed by the R4 strain is encoded by a well-studied ESBL 
gene, whose expression is not antibiotic inducible, one question 
remains to be addressed: how cell populations control CTX-M- 14 
gene expression, adapt to different antibiotic dosages, and thus 
display heteroresistance. If there is a regulator controlling hydro- 
lase gene expression, it must be a common cellular factor in all 
cells, since all hydrolase-possessing clinical isolates we have tested 
so far exhibited heteroresistance. Therefore, we ascribe to the 
growth rate the role of such a regulator, since it has been recog- 
nized as a global negative regulator of gene expression (26). When 
a cell dedicates itself to division and growth, generally speaking, 
gene expression other than cell reproduction is always repressed. 
On the contrary, when a cell stops dividing, transcription and 
translation become active and gene products accumulate. Mean- 
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FIG 5 Single-cell observation of phenotype transition under antibiotic stress, (a) Box charts showing average GFP intensities and growth rates of the R4 strain 
at the end of 2-h treatment in four escalating ceftriaxone concentrations. For each antibiotic dosage, 50 cells were randomly selected and observed, (b) The 
distribution of cellular GFP intensity versus the growth rate (frequency in color scale) of R4 cells after the same antibiotic treatment as shown in panel a. (c and 
d) Corresponding graphs for the S4 strain, (e) Dynamics of growth rate (GR) and GFP intensity (I) of two daughter cells that have departed from a single mother 
cell. The daughter cell CI shows an exponentially increasing GFP intensity and reduced growth rate, whereas the daughter cell C2 exhibits a relatively lower GFP 
intensity but a higher growth rate, (f) Scheme of growth-centric feedback loops that control the expression of hydrolase in a stringent way. 
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while, extensive expression of nonessential genes can also suppress 
growth due to an overconsumption of cellular energy and re- 
sources. Our results support the idea of cross-inhibition between 
growth and CTX-M-14 expression, i.e., all fast-growing cells ex- 
press smaller amounts of CTX-M-14 (borne as a nonessential 
gene by the cells), whereas all high-expression cells are slow grow- 
ing. There is a third subpopulation, where its cells are neither fast 
growing nor high expression (Fig. 5b). The reason behind this 
remains unclear; one of the possibilities is competitive genes that 
are highly expressed in the cells exerting resistance to antibiotic 
stress through mechanisms other than CTX-M-14. In addition, 
cells in the clonal populations vary their plasmid copy numbers 
not only from cell to cell but also from time to time. The growth 
rate dependence of expression is particularly pronounced for 
genes on a plasmid, which replicates independently of the chro- 
mosome and accumulates in high copy number in nondividing 
cells, leading to a collective gene dosage effect on expression. In 
fact, altering gene dosage is one of the major mechanisms bacteria 
used to regulate gene expression (26). This is confirmed by our 
results, in which plasmid copy number in the HIG cells is higher 
than that in the LIG cells. 

We are now able to put together negative feedback connecting 
three primary variables: antibiotic concentration, growth rate, 
and hydrolase expression (Fig. 5f). The net negative feedback cir- 
cuit helps the system to reach a balance; i.e., for each antibiotic 
concentration, there should be an adaptive growth rate and level 
of hydrolase expression. To test the regulation effect of the feed- 
back circuit on hydrolase gene expression and the growth rate, we 
compared the behavior of both the resistant (R4) and sensitive 
(S4) strains in their endurable ceftriaxone concentration ranges (0 
to 64 ng/ml for S4, corresponding to 0 to 2,048 ju-g/ml for R4) and 
found that S4 cells did not show an obvious gradient in both 
growth rate and GFP intensity (Fig. 5c); in other words, cells did 
not switch to the high-expression and slow-growing phenotype 
(Fig. 5d) without the link of hydrolase expression. 

Nested within the three-parameter circuit, growth rate and hy- 
drolase expression form double-negative feedback (Fig. 5f). In 
systems biology, such double-negative feedback conveys effects 
similar to those of positive feedback that augments gene expres- 
sion on little stimulation (27). In individual cells, expression of 
genes involved in positive feedback is triggered by intrinsic mo- 
lecular noise and leads to bi-/multistable expression (28). In our 
case, such a circuit may represent a major contributor to the quan- 
tized GFP expression observed in flow cytometry. A recent report 
demonstrated similar effects of growth-centric positive feedback 
that leads to heterogeneity in resistance to other categories of an- 
tibiotics and resistance mechanisms (29). Since the growth rate is 
a critical parameter sensitive to the intensity of environmental 
stresses, growth-centric control of gene expression may be a gen- 
eral part of the bacterial response to environmental stresses. Un- 
der some conditions, such as in biofilm, cells are found to be 
heterogeneous in growth rate, resistance, and extracellular poly- 
meric substance secretion (30), where growth control may play 
key roles as well in phenotypic heterogeneity and coordination of 
population-based strategies for adaptation and survival. 

In our study, the growth-centric feedback model provides a 
basic explanation for the bacterial population's response to envi- 
ronmental stresses. Through the adaptation circuit, the resistance 
gene product (CTX-M-14 in our case) is connected to the strength 
of environmental stresses and is accurately controlled without 



overconsuming cellular energy and other resources, thus helping 
the population as a whole to reach its maximal growth or better 
propagation. 

Bacterial cellular and population heterogeneity under environ- 
mental stresses has been observed and is proposed to be a result of 
genetic circuitry (31, 32), expression noise (33, 34), and asymme- 
try division (35, 36). Our current view suggests the central role of 
growth control, which can be fine-tuned by several cellular mech- 
anisms, including rRNA biosynthesis (37), plasmid copy number 
control (38), competitiveness in gene expression, and stringent 
control of cellular energy. Our study merely provides the simplest 
case of heteroresistance, where a single hydrolase gene is the only 
resistance mechanism. In the real world, however, clinical isolates 
may recruit multiple resistance genes or regulatory elements, such 
as those encoding outer membrane proteins or efflux pumps, 
which are essential genes of the host cells and are under stringent 
regulation as well (39). If a stronger positive feedback loop is pres- 
ent, expression of these resistant genes may remain constant un- 
der a wider range of antibiotic concentrations, even leading to 
on-and-off expression control (40). Thanks to the development of 
single-cell technology, such as cell-on-a-chip (23), single- 
molecule detection (41), and fluorescence time-lapse microscopy 
(42), real-time cell activity can now be traced precisely in both 
time and space, helping in deciphering the major controlling fac- 
tors of gene expression and growth. 

Our results provide a novel insight into mechanisms by which 
bacteria adapt to antibiotic stress and a path forward in under- 
standing antibiotic resistance assay and treatment. First, cellular 
and population heterogeneity of bacterial pathogens may interfere 
with clinical assays for drug resistance and provides a compelling 
explanation for poor repeatability of MIC or resistance judgment 
in some cases. Second, pathogens with weak hydrolase gene ex- 
pression may appear to be drug sensitive but exhibit strong resis- 
tance after antibiotic treatment that may occur rapidly and even 
lead to treatment failure. Application of resistance gene discovery 
techniques is of the essence to avoid false-negative results in clin- 
ical practice. 

Conclusions. In conclusion, our observations at single-cell 
resolution provide new insights into bacterial responses to antibi- 
otic stress. Bacterial progenies exhibit typical heterogeneity in 
both growth rate and stress-related gene expression and deal with 
environmental stresses in a collaborative way instead of being se- 
lected. Control of the growth rate provides opportunities for the 
expression of stress-related genes and is capable of serving as a 
critical regulator of these genes. Meanwhile, the growth rate pro- 
vides a primary response to environmental stresses and thus 
makes gene expression adaptive to stresses through stringent con- 
trol of the growth-centric feedback loops. These feedback loops 
provide a deeper understanding of how bacteria govern heteroge- 
neous phenotypes as an intrinsic strategy to reach maximal pop- 
ulation growth under variable environmental stresses. Therefore, 
we should always keep in mind that growth-controlled cellular 
and population heterogeneity may interfere with identification of 
resistant bacterial pathogens and antibiotic treatment in clinical 
practice. 

MATERIALS AND METHODS 

Clinical isolates, Etest, and PAP. All clinical isolates were collected in a 
multicenter epidemiological study for bacterial resistance in China from 
2006 to 2008. Etests (strips provided by bioMerieux) were used according 
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to the CLSI (Clinical and Laboratory Standards Institute) protocol and 
manufacturer's instructions with Mueller-Hinton agar plates (BD Difco, 
USA). The population analysis profile (PAP) assay was performed as pre- 
viously described (43). Briefly, a colony from an overnight plate of each 
strain was inoculated and cultured in fresh LB broth for 5 h, and then 
10 _1 , 10~ 2 , 10~ 4 , 10~ 5 , and 10~ 6 dilutions were prepared in saline. An 
equal volume (100 fil) of each diluted culture was spread onto LB plates 
with ceftriaxone or imipenem concentrations of a double-diluted range of 
0 to 1,024 (ng/ml. Colonies were counted after 48 h of culture, and the 
log 10 CFU/ml value was plotted against the antibiotic concentration. 
Ceftriaxone was purchased from Roche (People's Republic of China). 

Plasmid construction. The plasmid pUA66 (kindly provided by Uri 
Alon from the Weizmann Institute of Science in Israel) is constructed 
based on low-copy-number pSClOl and bears GFP-mut2 with a strong 
ribosome binding site as a reporter to measure the rate of transcription 
initiation (16). The W<Jctx-m-i4 genes and upstream sequences in eight 
clinical isolates were completely sequenced. They are 100% identical to 
that of the CTX-M-14 gene and its upstream sequence (GenBank acces- 
sion no. EU418916). The 276-bp upstream sequence of the CTX-M-14 
ORF was PCR amplified and inserted into BamHI of pUA66 as the pro- 
moter for GFP expression. The sequence of both the promoter and CTX- 
M-14 ORF was PCR amplified and inserted into the Bglll site. Then, the 
plasmid was introduced into DH5a competent cells (Invitrogen, Inc.), a 
common host cell for subclones derived from E. coli MG1655, and then 
cultured in LB medium. 

FCM and cell sorting. All flow cytometry (FCM) experiments were 
performed using a Becton-Dickinson FACSCalibur instrument (with ex- 
citation and emission light at 488 nm and 525 ± 12.5 nm, respectively, at 
a voltage of 650 V for GFP fluorescence). Wild-type DH5a was used as a 
GFP-negative control. About 10 5 to 10 6 cells (events) within a forward 
and side light scatter gate were assayed for the distribution of GFP expres- 
sion. FITC gates were used to sort cells with high or low GFP intensity. The 
sorted cells were collected in fresh LB medium for further analysis or 
culture. Original data were analyzed using the software program Flowjo 
(TreeStar, Inc.). 

Copy number determination of plasmid and CTX-M-14 transcript. 

DNA and RNA of the sorted cells were isolated using a commercial ex- 
traction kit (Qiagen bacteria; DNA minikit for DNA isolation and Am- 
bion PureLink RNA minikit for RNA isolation) . The relative copy number 
of the plasmid in each sample was determined by comparing the threshold 
cycle (C r ) value of the CTX-M-14 gene to that of the chromosomal gene 
gapA using quantitative real-time PCR. All amplifications were performed 
on an RT-Cycler amplifier (Capitalbio, China) with SYBR green (SYBR 
Premix Ex and TaqTM II; Takara). The relative mRNA of CTX-M-14 was 
determined using real-time PCR after being reverse transcribed to cDNA 
(PrimeScript RT reagent kit; Takara). Primers were as follows: for CTX- 
M-14, F,5' GGC TCA AAG GCA ATA CGA 3 ' ; R, 5 ' TTATCACCCACA 
GTC CAC GA 3 ' ; for gapA, F, 5' ACT TCG ACA AAT ATG CTG GC 3 ' ; R, 
5' CGG GAT GAT GTT CTG GGA A3'. 

Microfluidic chip and time-lapse microscopy. A polydimethylsilox- 
ane (PDMS) microfluidic cell-on-a-chip was fabricated with standard soft 
lithography, basically according to a previous design (44). Briefly, parallel 
thin rectangle culture chambers were used for observation of single-cell 
resolution. The chambers inside were supported by a pillar matrix and 
kept at a height of 1.5 u.m to constrain expanding E. coli colonies in 
monolayers. A feeding channel with a fast flow of fresh LB exchanged 
nutrition and wastes with the culture chambers by diffusion and con- 
trolled ceftriaxone concentrations as well. About 10 7 cells in exponential 
phase were resuspended in 50 /id fresh LB and introduced into the culture 
chambers. Images were obtained every 5 min at points of interests by 
using a Nikon Ti time-lapse microscopy instrument equipped with a PFS 
(perfect focus system) focus-maintaining unit and a charge-coupled- 
device (CCD) camera. The temperature of the device is controlled at 37 ± 
1°C. We processed time course images using the ImageJ software program 
and recorded the fluorescence intensity (AU) and cell length (/um). Be- 



cause the cells displayed a filamentous shape instead of dividing when 
treated with ceftriaxone, the doubling time could not be used to describe 
the growth rate in our study. We therefore defined the growth rate for each 
time point as dL/(dt ■ L), which means the elongation rate of each length 
unit (where d is differential and L is length of cell) . 
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